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The afc-plane thermal conductivity k of single-crystalline hexagonal MgB2 has been measured as 
a function of magnetic field H with orientations both parallel and perpendicular to the c-axis and at 
temperatures between 0.5 and 300 K. In the mixed state, k(H) measured at constant temperatures 
reveals features that are not typical for common type-II superconductors. The observed behavior 
may be associated with the field-induced reduction of two superconducting energy gaps, significantly 
different in magnitude. A nonlinear temperature dependence of the electronic thermal conductivity 
is observed in the field-induced normal state at low temperatures. This behavior is at variance with 
the law of Wiedemann and Franz, and suggests an unexpected instability of the electronic subsystem 
in the normal state at T w 1 K. 

PACS numbers: 74.7D.-b, 74.25.Fy, 74.25.-q 



I. INTRODUCTION 

The discovery of superconductivity in MgB2 below 
an unexpectedly high critical temperature T c of the or- 
der of 40 K 1 initiated intensive studies of its physical 
properties. Numerous results indicate that the super- 
conducting state of MgB2 is conventional in the sense 
that the electron pairing is mediated by the electron- 
phonon interaction. In most reports the superconduct- 
ing energy gap is claimed to be nodeless, compati- 
ble with s-wave pairing. However, various types of 
experiments, 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10,11 ' 12 ' 13 mainly using powder 
or polycrystalline samples and often surface-sensitive, 
have given evidence for two gaps of different magnitude 
in the quasiparticle excitation spectrum of this super- 
conductor. Calculations of the Fermi surface of this 
material 14 reveal three-dimensional sheets (7r-bands) and 
two-dimensional tubes (cr-bands) and it seems quite pos- 
sible that the gap values for these different parts of the 
Fermi surface differ substantially. It has been argued 15 
that the hole-like quasiparticles on the 2D parts of the 
Fermi-surface experience the larger superconducting gap 
with a maximum value close to 1.76fcsT c , as predicted by 
the original weak-coupling BCS theory. 16 The second and 
smaller gap is associated with the 3D sheets of Fermi sur- 
face. This intriguing situation and other possible anoma- 
lous features of this seemingly simple metallic compound 
ought to be checked experimentally on single crystalline 
material. 

Below we present the results of measurements of the 
thermal conductivity n and the electrical resistivity p 
parallel to the basal afr-plane of the hexagonal crystal 
lattice of MgB2 as a function of temperature T between 
0.5 and 300 K, and varying magnetic fields H between 
and 50 kOe, oriented both parallel and perpendicu- 
lar to the c-axis. In the mixed state, the observation of 
a rapid increase of the electronic thermal conductivity 
with increasing H at constant T <C T c is consistent with 
a field-induced reduction of the smaller of the two super- 
conducting energy gaps. Another important observation 
is that of an unusual nonlinear temperature dependence 



of the electronic thermal conductivity at T <C T c in the 
field-induced normal state, in obvious disagreement with 
the law of Wiedemann and Franz. 



II. SAMPLE AND EXPERIMENT 



The investigated single crystal with dimensions of 
0.5 x 0.17 x 0.035 mm 3 was grown with a high-pressure 
cubic anvil technique, similar to the one presented in 
Ref. 17 but with a slightly different thermal treatment, 
described in Ref. 18. The high quality of similar crystals 
from the same batch was confirmed by single crystal X- 
ray diffraction and energy dispersive (EDX) X-ray anal- 
ysis. A standard uniaxial heat flow method was used for 
the k(T, H) measurements. Temperatures between 0.5 
and 2.4 K were achieved in a non-commercial 3 He cryo- 
stat and the regime between 2 and 300 K was covered by 
using a 4 He flow cryostat. The thermometers for moni- 
toring the temperature difference in the k(T, H) measure- 
ments were Chromel-Au+0.07%Fe thermocouples with 
thin wires (25 /im diameter) for T > 2 K and a pair 
of ruthenium oxide thermometers below 2.4 K. In the 
region of temperature overlap, the results of the mea- 
surements of k(T) gave identical values for both types of 
thermometers. Since the thermopower of the Au+Fe al- 
loy is strongly field-dependent at low temperatures, spe- 
cial efforts were made to ensure a reliable calibration of 
the thermocouples in magnetic fields. Possible errors due 
to the thermal conductivity of connecting wires are es- 
timated to be below 1% of the total measured thermal 
conductivity. Additional measurements of the electrical 
resistivity p(T) in the afr-plane in magnetic fields oriented 
along the hexagonal c-axis were made as well. The elec- 
trical resistivity was measured by employing a 4-contact 
configuration with the same contacts used for measuring 
the voltage and the temperature difference, respectively. 
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FIG. 1: In-plane electrical resistivity p(T) of hexagonal 
MgB2 . The upper inset emphasizes the low temperature part 
for H = and 50 kOe oriented along the c-axis. The lower 
inset demonstrates the cubic temperature dependence of p(T) 
for T < 130 K. 




FIG. 2: Thermal conductivity vs temperature in the ab-planc 
of MgB2 in zero magnetic field. The solid lines represent k(T) 
measured on poly crystalline samples (data from: a - Ref. 25, 
b - Ref. 26, c - Ref. 23, d - Ref. 24). 



III. RESULTS 



A. Electrical resistivity 

The electrical resistivity p(T) in zero magnetic field is 
presented in Fig. 1. The narrow (AT C = 0.15 K) su- 
perconducting transition occurs at T c = 38.1 K. As it is 
shown in the lower inset of Fig. 1, in the temperature 
region between T c and 130 K, p(T) may well be approx- 
imated by p = po + AT 3 , where po and A are constants. 
At temperatures below 50 K, p ss po = 2.0 pfl cm. A 
cubic p(T) dependence has often been observed in multi- 
band transition metals and is associated with interband 
electron-phonon scattering. 19,20 The upper inset of Fig. 1 
emphasizes p(T) close to the superconducting transitions 
in fields of H = and H = 50 kOe, respectively. For our 
sample, the bulk upper critical field for the field direction 
parallel to the c-axis if£ 2 is about 30 kOe at zero tem- 
perature and decreases with increasing temperature, 21 
as shown in the inset of Fig. 3. Therefore the data for 
H = 50 kOe were obtained with the bulk of the sam- 
ple in the normal state. The abrupt slope change in 
p{T,H = 50 kOe) at 17 K is most likely related to the 
onset of spurious superconductivity in the surface region 
of the sample. 21,22 



B. Thermal conductivity 

The thermal conductivity data k(T) in zero mag- 
netic field are presented in Fig. 2. The k(T) values 
are about an order of magnitude higher than previously 
reported for polycrystalline samples. 23 ' 24,25,26 Also the 



overall temperature dependence of k is quite different 
from those earlier data. Instead of a monotonous in- 
crease with temperature we note a distinct maximum of 
k(T) at T ~ 65 K. The cause of these differences is obvi- 
ously the strong influence of intergrain boundaries on the 
heat transport in polycrystals, which masks the intrinsic 
mechanisms of quasiparticle scattering. As may be seen 
in Fig. 2, no anomaly in k(T) provides evidence for the 
superconducting transition at T c . A distinct change of 
slope in k(T) at approximately 6 K is observed, however. 

In Fig. 3 we display k(T) in the range between 0.5 and 
40 K, measured at selected fixed magnetic fields, oriented 
parallel to the c-axis. The practically overlapping curves 
for H = 33 and 50 kOe indicate that for this field ori- 
entation and these H values, the normal state has been 
reached in the whole covered temperature range. The ini- 
tial decrease of k with increasing field is relatively large 
at higher temperatures, but is progressively reduced with 
decreasing temperatures. It finally turns into an increase 
of k{H) at constant temperatures below approximately 
1 K. This behavior of the low-temperature thermal con- 
ductivity in the mixed state < H < H c2 is better il- 
lustrated in Fig. 4 where we present the k(H) curves, 
measured at constant temperatures below 8 K for field 
directions both parallel and perpendicular to the c-axis. 
The typical features of these curves are the rapid ini- 
tial decrease of k with increasing field, narrow minima 
in n{H) at field values that are low with respect to H c2 , 
and a subsequent s-shape type increase of n with fur- 
ther increasing field. The latter feature is particularly 
pronounced for the k(H) curves at the lowest tempera- 
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FIG. 3: Thermal conductivity vs temperature in the afo-planc 
of MgB2 for several values of magnetic fields parallel to the 
c-axis. The arrow indicates the zero-field critical temperature 
T c . The temperature dependences of the upper critical fields 
Hc2 and as established in Ref. 21 are shown in the inset. 



tures. The low-field increase is practically independent 
of the field direction as may clearly be seen in Fig. 4, 
where data for H \\ c (open symbols) and H _L c (solid 
symbols) arc displayed for comparison. The increasing 
slope at higher fields and H \\ c is certainly caused by 
approaching the normal state at H c2 . This trend is not 
observed for H _L c, for which H c2 is estimated to be 
about 130 kOe at low temperatures, obviously far be- 
yond our experimental possibilities. 



IV. DISCUSSION 

The thermal conduction of a superconductor is usually 
provided by electronic quasiparticles (n e ) and phonons 
(fv p h), such that 



(1) 



Upon decreasing the temperature to below T c in zero 
magnetic field, the reduction of the number of unpaired 
electrons leads to a decrease of n e and an increasing K p h- 
The overall behavior of n(T) in the superconducting state 
depends on the relative magnitudes of K e and K p h and 
also on the strength of the electron-phonon interaction. 
Applying external magnetic fields induces vortices in the 
sample. The quasiparticles associated with the vortices 
not only provide additional contributions to phonon scat- 
tering and hence a reduction of K p h, but also enhance K e . 
The competition of these two processes leads to k(H) 
curves as shown in Fig. 4. 
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FIG. 4: Thermal conductivity in the basal plane of MgB2 vs 
H at several fixed temperatures. The arrows denote the upper 
critical field H C 2 for H \\ c. The closed and open symbols 
correspond to the field direction perpendicular and parallel 
to the c— axis, respectively. 



In what follows, we first discuss the low-temperature 
temperature part of n{T) in the field-induced normal 
state, i. e., for H > H c2 . Subsequently we discuss n{T) 
in the superconducting state for H = and finally turn 
to the data of k(H), measured at constant temperatures 
in the mixed state, i. e., H < H c2 . 



A. Thermal conductivity in the normal state 

One of the main problems encountered in the analy- 
sis of the thermal conductivity of a conductor is that a 
separate identification of the two terms in Eq. (1) at arbi- 
trary temperature is not straightforward. In the normal 
state, a convenient and often used way to estimate K e 
is to employ the Wiedemann- Franz law (WFL), relating 
the electrical resistivity and the electronic contribution 
to the thermal conductivity via 



n e (T) - L T/p(T), 



(2) 



where L = 2.45 x 10~ 8 W £1 K~ 2 is the Lorenz num- 
ber. The validity of this law requires an elastic scatter- 
ing of electrons and it is well established that Eq. (2) is 
applicable if the scattering of electrons by defects dom- 
inates. This is usually true at low temperatures, where 
p(T) w p . In our case, the data in Fig. 1 suggest 
that Eq. (2) is applicable at temperatures below about 
50 K. Inserting the appropriate numbers into Eq. (2) sug- 
gests that, at temperatures T c < T < 50 K, K e (T) pro- 
vides about half of the total thermal conductivity. As 
we demonstrate below, the applicability of the WFL for 
MgB 2 at low temperatures is questionable, however. 

Before discussing the validity of WFL at low tempera- 
tures, we note the complications that arise from unusual 
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features in the p{T) data shown in Fig. 1, most likely 
caused by the influence of superconductivity in minor re- 
gions of the sample with enhanced T c and H C 2, different 
from the corresponding bulk values. The onset of this su- 
perconductivity is clearly seen in the upper inset of Fig. 1 
by inspecting p(T) measured in a field of 50 kOe, substan- 
tially exceeding i? c2 (0), the maximum bulk upper critical 
field of the same sample. 21 Superconducting traces in a 
minute fraction of the sample may cause a considerable 
reduction of the total measured electrical resistivity but 
leave the thermal conductivity virtually unchanged. In 
this case a failure of the WFL would not be surprising. 
However, because the electrical resistivity of the sam- 
ple above T c is practically temperature- independent, wc 
can, with a great deal of certainty, expect that the intrin- 
sic electrical resistivity of the bulk remains constant also 
at lower temperatures. The constant residual resistivity 
po is caused by defects which actually set the maximum 
mean free path for electrons, independent of temperature. 
Indeed, results by Xu et al. 27 of p(T) measurements on 
single crystalline MgB 2 at magnetic fields which presum- 
ably exceed the upper critical field of the minor phase, 
demonstrate that in, e.g., H = 90 kOe the electrical re- 
sistivity remains constant with decreasing temperature 
down to at least 2 K. 

In Fig. 5, we plot as a solid line, the normal-state elec- 
tronic thermal conductivity k^ fl below 8 K, calculated 
using Eq. (2) with the experimental value po = 2.1 pflcm 
for H = 33 kOc. The measured total thermal conductiv- 
ity, shown by open circles, is considerably higher than 
k^ fl across the entire covered temperature range. For 
the estimate of the upper limit of the phonon contribu- 
tion we assume that the minima of k(H) shown in Fig. 4 
are caused by the competition of a decreasing K ph and an 
increasing K e . With this interpretation it is clear that the 
values of K m i n (H) represent at most the maximum value 
of the lattice contribution K p h- The smooth interpolation 
between these minimum values of k{H) for different tem- 
peratures, denoted as K^jf*, is shown as the broken line 
in Fig. 5. The difference between the measured thermal 
conductivity at H = 33 kOe and the upper limit of the 
phonon contribution, nf m = n — k™£ x , obviously repre- 
senting the lower limit of the electronic contribution, is 
shown in Fig. 5 by open squares. It may be seen that 
at least below 8 K, the electronic contribution is consid- 
erably larger than predicted by WFL. With increasing 
temperature, K™ m approaches the WFL prediction, as is 
demonstrated in the inset of Fig. 5. Nevertheless, it is 
impossible to identify the temperature limit where the 
validity of the WFL is recovered if such a limit exists 
at all. Considering our procedure it is clear that the 
true electronic contribution exceeds our estimate, partic- 
ularly towards the upper end of the considered temper- 
ature range. However, at very low temperatures where 
k(H = 0) <C K n (e. g., at 0.60 K the zero field thermal 
conductivity is less than 6% of the normal-state thermal 
conductivity), nf in must be very close to the true K e . 
Thus the observed deviation from the WFL and in par- 



ticular its temperature dependence, revealing a peak-like 
structure of n e /T vs T, shown in the inset of Fig. 5, is a 
reliable result of our investigation. 

The violation of the Wiedemann- Franz law at low tem- 
peratures is very unusual because the validity of this law 
is expected to hold for the Fermi-liquid ground state of 
common metals. A similar set of data of Hill an cowork- 
ers for (Pr,Ce)2Cu04 was interpreted as evidence for 
a breakdown of the Fermi-liquid theory for this oxide 
material. 28 The non-Fcrmi-liquid behavior of cupratcs 
was speculated to be the consequence of possible spin- 
charge separation, a scenario that is considered in rela- 
tion with high-T c superconductivity. 29 The same argu- 
ments are hardly relevant for the case of MgB 2 where a 
spin-charge separation is not expected. Another expla- 
nation for the results of Ref. 28 was recently offered in 
Ref. 30. It was argued that a peak-type structure of n e /T 
plotted versus T might be the consequence of some kind 
of transition leading to a gap formation in the electronic 
excitation spectrum of the normal state. This scenario 
would lead to a peak of n e /T at T , close to the transition 
temperature, an exponential decrease of n e /T below T , 
and a gradual approach to the Wiedemann-Franz value 
of L a /p at T > To. Such a peak at T ~ 1 K indeed 
follows from our analysis of the MgB 2 data (see the inset 
of Fig. 5). Unfortunately our data set does not extend 
to low enough temperatures in order to confirm the ex- 
ponential temperature dependence of n e well below T . 



B. Thermal conductivity in the superconducting 
state (H = 0) 

A rather unexpected feature in the temperature depen- 
dence of k in zero magnetic field is the absence of even the 
slightest manifestation of the transition at T c — 38.1 K, 
as has already been mentioned in previous reports on 
k(T) for polycrystallinc materials. 23 ' 24,25,26 This obser- 
vation is quite unusual for superconductors with non- 
negligible phonon heat transport, because the opening 
of the superconducting gap rapidly reduces the rate of 
phonon scattering on electrons and should lead to a fast 
increase of « p h below T c . The assumptions that either the 
phonon-electron scattering is much weaker than phonon- 
defect scattering, or that « p h is negligibly small in the 
vicinity of T c , which, in principle, might explain the ab- 
sence of a k(T) feature at T c , are all incompatible with 
the observation that applying a relatively weak external 
magnetic field of 0.63 kOe, introducing some additional 
quasiparticles in the cores of vortices, considerably re- 
duces the thermal conductivity at intermediate tempera- 
tures (see Fig. 3). The possibility that the enhancement 
of Kph below T c is exactly compensated by a reduction 
of K e , is not considered because, for a BCS supercon- 
ductor, the slope change in n p h{T) at T c is much more 
pronounced than that in K e (T). 31 The latter statement 
is not true for special cases of extremely clean samples of 
strong-coupling superconductors, such as Pb and Hg, 32,33 
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FIG. 5: Normal-state thermal conductivity measured in 
H = 33 kOe between 0.6 and 8 K (open circles). The upper 
limit of the phonon thermal conductivity K™ h ax (dashed line) 
and the lower limit of the electronic contribution ft™ ln (open 
squares) are calculated as described in the text. The solid 
line represents the electronic contribution k^ vfl calculated 
using the Wiedemann-Franz law in Eq. (2). In the inset, we 
plot K™ ln (open squares) and Z/o/po (solid line), respectively. 
The error bars in the inset mark the maximum uncertainty of 



where the scattering of electrons by defects is negligibly 
small. Our sample of MgB 2 with p(300)/p(0) « 6.8 can- 
not really be regarded as fulfilling the extreme clean limit 
criteria. It is possible, however, to account for both the 
absence of a feature in k(T) at T c and the obvious slope 
change in k(T) centered around 6 K, by postulating that 
the superconducting energy gap A(T) for quasiparticles 
which strongly interact with low- frequency phonons, is 
considerably smaller than the values given by the BCS 
theory. 

In the simplest approximation, the phonon thermal 
conductivity can be calculated as 



K p h = (« 2 /3) J C(u))t(u))(1u), 



(3) 



where v is the mean sound velocity, and C(u>) and t(u>) 
are the specific heat and the average relaxation time of 
a phonon mode with frequency uj, respectively. The to- 
tal phonon relaxation rate may be calculated by assum- 
ing that the simultaneous influence of all independent 
phonon scattering mechanisms n can be accumulated in 
the form 



" 1 = E- 



The phonon-electron relaxation time r p _ e changes most 
rapidly upon the opening of the superconducting gap. 



In Ref. 31, the phonon-electron scattering time in the 
superconducting state takes the form 



r p s _ e H-.g(c/T,A/T)r;_ e H, 



(5) 



where T™_ e (aj) is the normal-state relaxation time. The 
function g(u>/T, A/T) is quite complicated, but its main 
feature is a step-like increase of the phonon relaxation 
rate at the phonon frequency lo — 2A/K, a consequence 
of the fact that a phonon with an energy less than 2A 
cannot brake a Cooper pair and interacts only with quasi- 
particles that have already been excited above the gap. 34 
In the "dominant phonon" approximation 35 it is assumed 
that, at temperature T, the main contribution to the heat 
transport in the lattice is due to phonons with frequen- 
cies close to Wdom, where ?iWdom ~ 3.8fcsT. For weak- 
coupling BCS superconductors, where A(0) = 1.76fcsT c 
and A(T) is a standard function tabulated in Ref. 36, 
hu>d om (T) = 2A(T) is always fulfilled at the temperature 
0.73T C , i. e., not far below T c . This is the reason why 
K P h(T) increases rapidly below T c , leading to a maximum 
of the measured k, typically close to T c /2. 34 For MgB2, 
instead of a peak near T c /2, k(T) exhibits a distinct fea- 
ture at about T c /6 ~6K, which can also be regarded as 
a peak-type structure on top of the background which de- 
creases monotonously with decreasing T. This suggests 
that the relevant superconducting energy gap is equal to 
TiWdom at much lower temperatures than 0.73T C . Hence, 
for the quasiparticles which scatter phonons most effec- 
tively, the energy gap A(T) is about 3 times smaller than 
the values given by the original BCS theory. 

In principle, more information could have been ex- 
tracted from our data by direct comparison with ex- 
isting theories for the thermal conductivity in multi- 
band superconductors. 34 ' 37,38 ' 39 However, since the WFL 
seems to be invalid in the normal state of MgB2 at low 
temperatures, any attempt to analyze k(T) quantita- 
tively in the superconducting state is hampered by the 
difficulties in separating n e and n p h in a reliable manner. 



C. Thermal conductivity in the mixed state 

Before discussing the features of k in magnetic fields, it 
is important to note that the zero-field values of k at tem- 
peratures T <C T c are almost entirely due to the phonon 
contribution. The electronic thermal conductivity in the 
superconducting state K e>s {T) can be estimated using the 
theory of Bardeen, Rickayzen, and Tewordt 31 . In their 
model 



where 



W f(y) 



K e,s — K e,nf(lj)j 



2F 1 (-y) + 2yln(l + e-y) + 



(6) 



2*1(0) 



and 

m 



-, (T<T C ) 
(7) 



1, (T > T c ), 
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as well as F n (-y) = J* °° z n (l + e z+v )- x dz, y = 
A(T)/kBT. Using for n e , n the values of shown 
in Fig. 5, and taking into account the lowest previously 
claimed value of A(0) = 1.7 meV, 10 it may be shown 
that K e ,s is negligibly small below 4 K, and therefore, 
k(H = 0) « Kp h - 

As it may clearly be seen in Fig. 4, the n{H) curves re- 
veal very similar features at all temperatures below 8 K. 
The most intriguing aspect of these curves is the very 
rapid increase of k at relatively weak fields, after the ini- 
tial decrease of K p h- This increase is undoubtedly due to 
a field induced enhancement of the number of electronic 
quasiparticles and thus to an increase of K e . For com- 
mon type II superconductors in their mixed state, the 
features of K e (H) are expected to depend on the ratio 
between the electron mean free path I and the coher- 
ence length £o- A rough estimate of t can be obtained 
from the residual resistivity po employing the Drude re- 
lation po = 3/[N £vFe 2 }, where No is the electronic den- 
sity of states at the Fermi level and vp/3 is the aver- 
age in-plane component of the Fermi velocity. Using 
the values N — 0.7 states/(eV unit cell) (Ref. 14) and 
vp = 4.9 x 10 7 cm/sec (Ref. 4), we obtain £ « 80 nm. 
Since this value is considerably larger than the in-plane 
coherence length £ a ^ = 11.8 nm (Ref. 21), our sample 
is in the moderately clean limit. In this limit, n e is ex- 
pected to be small at all fields below H c2 , except close to 
H c2 , where it grows according to 

Ke = « e , n (l-C T (ff C 2-H) 1/2 ), (8) 

where n oc T is the normal-state electronic ther- 
mal conductivity above H C 2 and Ct is a temperature- 
dependent coefficient. 40 As may be seen in Fig. 4, for 
MgB 2 , the largest positive slope dn/dH is observed well 
below H C 2 for both orientations of the magnetic field. Al- 
though the form of Eq.(8) implies that k(H) should scale 
with the value of LT C 2, the thermal conductivity depends 
only weakly on the field direction up to approximately 
6 kOe. This is amazing because the anisotropy of the 
upper critical field H^/H^ ~ 4.2 at low temperatures. 21 

For H _L c, after a steep initial increase, k(H) reaches 
a region where it exhibits only a relatively weak in- 
dependence. The same tendency is also observed for H\ |c 
but it is partly masked by yet another increase of k(H) 
close to H c2 . For H ± c and T < 8 K, H$ ~ 130 kOe 
(Ref. 21 ) and is not accessible in our experimental setup. 
Therefore the region of weak ii-dependence extends to 
the highest fields reached in this study. 

The field dependence of K e , although very different 
from what one would expect from a conventional super- 
conductor, can qualitatively be explained in terms of a 
two-band model with two energy gaps of different magni- 
tude associated with each band. Nakai and co-authors 41 
analyzed such a model where one band with strong pair- 
ing (L-band) is responsible for superconductivity, and su- 
perconductivity in the second band (S-band) is induced 
by Cooper pair tunneling. Consequently, the two bands 
are characterized by a smaller gap As and a larger gap 



Al, and normal-state electronic densities of states at the 
Fermi level A^,s and Nq^, respectively. The analysis 
presented in Ref. 41 shows that the quasiparticle states 
in the vortices are highly confined in the L-band but only 
loosely bound in the S-band. Therefore the quasiparti- 
cle states of the vortices in the S'-band start to overlap 
already in weak fields and the resulting density of states 
equals that of the normal-state A*o,s at H -C H C 2- The 
situation can be visualized as a vortex lattice involving 
the L-band states, coexisting with the normal state in 
the S'-band where the energy gap is suppressed. This 
model explains very well the behavior of the electronic 
specific heat in a magnetic field. 2,42 The field-induced 
suppression of the smaller gap is claimed to be consistent 
with the results of point-contact measurements 3 ' 9 and re- 
cent scanning tunneling spectroscopy experiments. 43 In 
terms of the two-gap model, the saturation of the ther- 
mal conductivity much below H c2 may be regarded as 
the result of the closing of the energy gap in the S-band. 
The heat transport via quasiparticles of the band associ- 
ated with the larger gap is significant only in the vicin- 
ity of H c2 , and, eventually, above H c2 the full normal 
state electronic heat transport is restored. The lack of 
a substantial dependence of K e on the field orientation 
for H <C H C 2 is an obvious consequence of the weakly 
anisotropic 3D nature of the 7r-bands. From this we con- 
clude that the smaller gap must open in the 7r-band. The 
rapid increase of the number of quasiparticles in the tt- 
bands also naturally explains the very fast drop of K p h(H) 
in small fields, because the corresponding excited quasi- 
particles are the dominant scattering centers for phonons 
at low temperatures. 

A more quantitative analysis of n e {H) can be made 
for the lowest temperatures where, as may be seen from 
Fig. 3, the phonon contribution is relatively small in 
comparison with the field-induced electronic contribu- 
tion. At very low temperatures, the phonon scattering 
by electronic quasiparticles is less effective, 35 therefore 
K p h{H > 0) should not much deviate from K p h(H = 0). 
Indeed, at temperatures of 0.60 and 1.02 K, there is no 
initial decrease of n(H) in small fields. Assuming that the 
phonon contribution is essentially LLindependent and 
that the smaller gap is completely suppressed in fields 
exceeding 20 kOe, we may establish the individual con- 
tributions to K e of the quasiparticles associated with ei- 
ther the a- or the 7r-band. This is illustrated in Fig. 6. A 
possible reduction of K p h with increasing H could slightly 
change this ratio in favor of n e ^, but only by a few per- 
cent. The ratio n e ^jn e , a is 0.57/0.43, as estimated from 
the k{H) data at T — 0.60 K. This ratio is remarkably 
close to the ratio of the densities of electronic states in the 
two bands N 0i7r /N , a of 0.58/0.42, as calculated by Liu et 
al. 15 and 0.55/0.45 by Belashchenko et al. 44 . Similar ra- 
tios of 0.55/0.45 and 0.50/0.50 have been extracted from 
tunneling spectroscopy measurements 43 and from specific 
heat experiments, 2 respectively. From this comparison, 
we reach the important conclusion that the electron mean 
free paths on different sheets of the Fermi surface are 
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close to being equal. Indeed, the electronic thermal con- 
ductivity may be calculated from n e = C e vp£/3, where 
C e is the electronic specific heat. Since C e> i oc Noj, and 
the afr-components of the Fermi velocity Vf,% (i — ft, o) 
are similar for different sheets of the Fermi surface in 
MgB 2 , the equality K e ,n / ' K e ,a ~ Nq^/Nq^ is tanta- 
mount to saying that, at low temperatures, the ratio of 
the electron mean free paths in different bands £ v /£ a is 
close to unity. This observation is essential in view of 
the current discussion of the possibly different impurity 
scattering rates in different bands of electronic states of 

MgB2 46,47,48 

Although the absence of any particular feature at T c in 
the zero-field n(T) data, discussed in Section IV B, gives 
only qualitative support for the existence of parts of the 
Fermi surface with a gap much smaller than predicted by 
the standard BCS theory, the magnetic field induced vari- 
ation of the low-temperature thermal conductivity may 
be regarded as strong evidence in favor of the multigap 
scenario. 

At the same time we believe that the model of Haas and 
Maki 49 for explaining thermodynamic and optical prop- 
erties of MgB2 is not appropriate. They proposed the 
k-dependence of a single energy gap to adopt the form 
of a prolate ellipsoid. Our claim is based on the com- 
parison of our data set with similar results for materials 
with strongly anisotropic gap functions. In Fig. 7, we 
redraw a figure from Ref. 50, which compares the field- 
induced variation of n e at temperatures well below T c for 
different conventional and unconventional superconduc- 
tors, amended by our data for MgB 2 at T = 0.60 K. The 
Nb data 51 reveal the typical salient features of a clean, al- 
most isotropic s-wave superconductor, and confirms the 
very weak energy transport by quasiparticles far below 



FIG. 7: The electronic thermal conductivity normalized to 
its normal state value vs. H/H C 2- The data for MgB2 are 
from this work, the results for Nb, UPt3, and LUN12B2C are 
from Refs. 51, 52, and 50, respectively. 



H c2 , in agreement with Eq. (8). A considerably faster, 
almost linear in iJ, increase of n e is observed for a su- 
perconductor with nodes in A(k), here exemplified by 
UPt3. 52 A similar H- variation of K e has been observed 
for LuNi 2 B 2 C, which led the authors of Ref. 50 to claim 
an anisotropy of the energy gap A max /A m j n > 10. The 
increase of n e (H) in MgB 2 is much faster than in any of 
these materials. For the field direction H _L c, more than 
half of the normal-state thermal conductivity is restored 
already at H = 0.05H c2 \ This means that, with increas- 
ing H, instead of the gradual increase of the number of 
quasiparticles contributing to the heat transport, that is 
characteristic of single-gap anisotropic superconductors, 
we are dealing with the suppression of an energy gap on 
a significant portion of the total Fermi surface by rela- 
tively weak magnetic fields. Thus the features of n e (H) 
of MgB 2 displayed in Fig. 7 may be regarded as a natural 
consequence of the existence of two different gaps. 



V. SUMMARY 

Our k(T, H) data provide evidence for a rapid field- 
induced enhancement of quasiparticles in the supercon- 
ducting state of MgB 2 well below H c2 , consistent with an 
efficient field-induced closing of the smaller energy gap, 
thus provoking a fast growth of the electronic thermal 
conductivity. At higher fields, the growth of K e (H) tends 
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to saturate until, in the vicinity of H C 2, the contribu- 
tion to k from the electrons of the band associated with 
the larger energy gap grows rapidly, merging into the 
practically field-independent thermal conductivity in the 
normal state above H C 2- 

At low temperatures, the electronic thermal conduc- 
tivity of the field-induced normal state is nonlinear in 
T and deviates considerably from the prediction of the 
Wiedemann-Franz law. This deviation peaks at about 
1 K, suggesting the existence of some transition provok- 
ing a gap formation in the electronic excitation spectrum 



close to this temperature. 
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